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Abstract 

A new, exactly solvable, Barbieri-Remiddi like equation for bound states of two scalar 
constituents interacting with massless vector particles is presented, both for stable and 
unstable particles. With the help of this equation the bound state spectrum is calculated 
to O(a^) for a SU(N) nonabelian gauge theory. The result for the abelian case reproduces 
the known result from the Foldy-Wouthuysen calculation. It is shown how different graphs 
as in the fermionic theory contribute to the spectrum to this order. Furthermore the bound 
state correction to the decay width for a weakly decaying system is calculated. This result 
is equal to its fermionic counterpart. Thus the theorem on bound state corrections for 
weakly decaying particles, formulated previously for fermions only, has been extended to 
the scalar theory. 
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1 Introduction 



While the discussion of fermionic bound states has a long history ||l|, much less attention 
has been paid to the similar problem with scalar constituents. Only the ladder approxi- 
mation with scalar interaction is a well known example and has already been discussed in 
the 50-s and 60-s in the framework of the Bethe-Salpeter equation. Indeed, to this day 
the only known fundamental matter fields are fermionic. But in supersymmetric theories 
for each fermion two scalar partners are required. Since some of them, probably stop or 
sbottom, could have masses within the reach of the next generation of e~^e~ accelerators, 
even the observation of bound states of those particles seems possible. These objects and 
systems built of scalar composite particles in atomic physics underline the need of an 
equally clear and transparent approach as the one developed for the fermionic case |Q. 
A recent attempt in this direction 0] splits the boson propagator in a particle and anti- 
particle propagator in order to be able to treat them like fermions. The spectrum is then 
obtained by constructing the Hamiltonian via a Foldy - Wouthuysen transformation and 
a perturbation theory a la Salpeter. This approach does not show significat advantage 
over the pure Foldy- Wouthuysen approach |Q and suffers also from the drawback that it 
will break down in higher orders due to the appearance of higher powers in the spatial 
momentum p. The Coulomb field appears in this formalism as an external field which 
makes this formalism not very reliable-looking. All these drawbacks can be circumvented 
by developing an exactly solvable zero order equation and subsequently using a systematic 
perturbation theory. 

To the best of our knowledge there exists no attempt in the literature to construct a 
solvable zero order equation for the BS equation containing two charged scalars interacting 
via a vector field. This goal will be achieved in section 2. 

In section 3 we will review briefly the BS perturbation theory and use it to calculate 
the spectrum of bound states for scalar particles with equal mass, both for the abelian and 
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nonabelian case to 0{a'^). This will be of importance if the stop has a narrow width. If the 
width becomes comparable to the level splittings this considerations can be understood 
as a determination of the scalar-antiscalar potential. 

The decay width is also subject of the second application we present in section 4. We 
calculate the bound state correction to the decay width T of system of scalar constituents 
to 0{a^T). 

Finally section 5 is devoted to the conclusions and to the discussion of our results. 

2 A bound state equation for scalar particles 
2.1 Stable particles 

As starting point we present here an exactly solvable equation for stable scalar particles 
which interact via a vector field. 

We start from the BS equation for a bound state wave function x 

where 5* is the exact scalar propagator, and K is the sum of all two scalar irreducible 
graphs. Both are normalized to be Feynman amplitudes. Furthermore, we have introduced 
relative momenta p and p', a total momentum P = Pi — P2, and we choose the center of 
mass (CM) frame where P = (Pq, 0) = (2m + E,6). 

As a first approximation to eq. (0) we would like to use beside the free relativistic 
scalar propagators the kernel due to the Coulomb interaction 

Kc{p,p)=ATTa — — — . 2 

[p-p'r 

For a nonablelian theory with gauge group SU{N) we use 



In this case x has to be a singlett in order that Kc represents an attractive force. The 
kernel @) has the drawback that it is po dependent and the exact solution of eq. with 
(0) is not known. However, in the nonrelativistic regime by the scaling argument [^] 

Pq ~ 0(ma;^), \p\ ~ 0(mQ;), (4) 
Po ~ 2m-0{ma^) 

we can start with an instantaneous approximation to the kernel since Pq is of O(a^m) in 
this region and may be included in the corrections afterwards. Doing this, we can perform 



the zero component integration on the propagator ( Ep = yw?+p^) 

dpo 1 



—I 



2vr [(f + Po)' - El + ze] [(-f + p,y - E^ + le] 
1 



2Ep - Po 2Ep + Po 
and it is quite easy to show that 



Epi^El - PS) 



(5) 



Ko{p,p')=47ra V (6) 
q ^ 

gives a solvable equation with the normalized solutions 

/^(P2 - 4^2) 

""^^^ " 'x/2i^[(f +Po)^ - El + .e][(-f +po)^ - El + 

= -e) (8) 

to the eigenvalues 

Po = Mf^ = 2m^l - a„ = £. (9) 

Eq. (|) is dictated by the requirement that x should acquire the same analytic properties 
as the underlying field correlators 

Xip) = /e^^^(0|T$t(^)$(-|)|P„), (10) 
Xip) = /e-^^^(P„|T$(^)$t(_|)|o>. (11) 



Using the integral representation for the step function which is included in the time 
ordered product, one derives eq. (H). 

Taking the equation for the Green function 



tGo = -Do + DoKoGo, (12) 

with 

r. ^ {2n)'S\p-p') 

instead of that for the BS wave function and using again (0) we find 



with 



and 



G„ = -F(p) ^"'ff-^^V (p-) (14) 



E^^^ (15) 



F(n) = V " ^ (IQ) 

Gc denotes the well known Coulomb Green function in momentum space. These solutions 
can be used for a systematic BS perturbation theory for scalar constituents, as will be 
demonstrated in the next section. 



2.2 Unstable Particles 

As has been shown recently by the author []^ for the fermionic case, an important sim- 
plification can be achieved in some calculations if the width of the bound state is already 
included in the zero order equation. Furthermore, if the width becomes comparable to 
the level shifts, this approach even becomes indispensible. For the scalar case this can be 
done by the replacement 

Ep ^E^ - tTm. (17) 
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While (p!7D leads to expressions for the BS wave functions which contain unpleasant ex- 
pressions for the particle poles it has the advantage that the propagator has the form as 
expected from the phase space of an unstable particle. Furthermore the above calculation 
remains essentially unchanged if we define the square root in ( p^Tf ) to be that with the 
negative imaginary part (clearly we demand F > and m > 0). Only the energy in the 



resulting equation for the Green function and thus in (14) changes to 



f, P2 _ 4^2 

4m 



The eigenvalues for Pq are 



Pon = 2mW 1 - a2 - z— ^ 2m - ma^ ^ + iV - 19 

\ m 4 4m 2 

In the case of the fermions we managed to construct wave functions independent of 

F. This was possible because the small components of the propagator containing Pq — iV 

instead of Pq + were projected away by the choice of an appropriate kernel K. This 

cannot be achieved in the scalar case and thus, surprisingly enough, the scalar wave 

functions look more complicated than the fermionic ones. A version for a zero order 

equation for decaying particles where the propagator is chosen in close analogy to the 

fermionic case has been developed in . In our present work we, instead, proceed in the 

spirit of our generalized approach. 



3 Perturbation Theory 

Perturbation theory for the BS equation starts from the BR equation for the Green 
function Go (eq- (HID) of the scattering of two fermions [^j which is exactly solvable. is 
the product of two zero order propagators, Kq the corresponding kernel. The exact Green 
function G may be represented as 

1 oo 

G = E xI!t, ^X^f + Greg = Go j:iHGor, (20) 
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where the corrections are contained in the insertions H and Greg is the part of G regular 
at Pq = Pn- It is easy to show that H can be expressed by the full kernel K and the full 
propagators D: 

H = ~K + Ko + iD-^ -iD^K (21) 

Thus the perturbation kernel is essentially the negative difference of the exact BS-kernel 
and of the zero order approximation. 

Expanding both sides of equation (^) in powers of Pq — Pn, the mass shift is obtained 

EM- 

AT 

AM-z— = {ho){l + (h,)) + {hog.ho) + Oih^). (22) 
Here the BS-expectation values are defined as e.g. 

^^^^^ ^ J J '^^^^^^P^^^^'^^'^P^P">^^'^'^P^^ ^2^) 
We emphasize the four-dimensional p-integrations which correspond to the generic case, 
rather than the usual three dimensional ones in a completely nonrelativistic expansion. 
We distinguish these two case by introducing the notation ((...)) for a four- dimensional 
expectation value and (...) for the usual nonrelativistic expectation value 

Of course, ( 123|) reduces to an ordinary "expectation value" involving d^p and $(p), when- 
ever h does not depend on po and pg- 

In (|22|) hi and Qi represent the expansion coefficients of H and Go near the pole at P„, 
respectively, i.e. 

oo 

H = J2hM-Pnr (25) 

m=0 

oo 

Go = Y.9m{Po-Pnr-' (26) 

m=0 

Similar corrections arise for the wave functions P, |10 



X^'^ = {9iho + l{h))x^'^ (27) 
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3.1 Fine structure 

As an application of this perturbation theory as well as of the new zero order equation 
for scalar particles developed in the last section, we will present here the calculation of 
the fine structure of two stable scalar particles interacting via a vector particle. Existing 
calculations |Q rely on a mix of Fouldy-Wouthuysen transformation and the iterated 
Salpeter perturbation theory. Our present approach is much more transparent and allows 
in principle the inclusion of any higher order effect in a straightforward manner. First we 
will calculate the fine structure for two scalars of equal mass interacting by an abelian 
vector field. Then we consider also the nonabelian case which could be of interest for the 
stop-antistop system. In this case we will calculate the spectrum up to order aj. 

Since in the zero order equation we have replaced the exact one Coulomb exchange 
(0) by Kq as given in we have now to calculate the contribution of —Kq + Kq to the 
energy levels. This is shown in fig. (pa . With 

ii-^c)) = -4-«y (^(^x(i>) x(p) = 

= -( ° p^^ j^Z^) = (28) 
, / 2m 0",? \ Ana , 



am = ^^{i^) (29) 

we obtain 

AMc := i(-^c + Ko)) = f{^) = ^^. (30) 

The fact that the p- integrations are well behaved and the result is of 0{a^) proves the 
usefulness of our zero order kernel. 

The transverse gluon of fig. |l]b gives rise to a kernel 

HT = -^[{p + Pj ^ 1 (31) 




a 



b 




WVA; 



transversal gluon 
Coulomb gluon 
transversal or Coulomb 



c 



Fig. 
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Performing the zero component integrations exactly and expanding in terms of the spatial 
momenta one obtains to leading order (c.f. HTlH ) 



Due to the fact that scalars can only form spin zero bound states, the the annihilation 
graph into one gauge particle (with spin one) contributes only for p-waves and thus is 
supressed by two additional powers in a. Furthermore, as in the fermionic case, it vanishes 
for the nonabelian theory due to the color trace since the bound states are color singlets. 
As can be seen from the above results the contribution of the transverse gauge field 
is equal for fermions and bosons. However, the relativistic correction to the Coulomb 
exchange appears to be different. Let us therfore check the contribution of this Coulomb 
correction from second order perturbation theory (fig. |^.a). These contributions give only 
rise to O(a^lna) effects in the fermionic theory. Since the leading Coulomb singularity 
is cancelled we may hope that we can replace the Green function by the free propagator. 




(32) 



(33) 
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Indeed it can be shown that the next terms of the Green function give only higher order 
contributions. 

Due to the presence of the zero component momentum in the scalar- Coulomb gluon 
vertex we observer that the contribution from fig.|[b diverges linearly. However, it is an 
easy exercise to show that in the sum of graphs (fig.|^.b + fig. this linear divergence 
cancels. Thus we regularise all the single graphs, sum up and find a finite result. We have 
used dimensional regularization as well as a one dimensional Pauli-Villars regularization. 
Both give the same result for the finte parts of the intgrals. 

AM,,, = {{h^^^)) + {{{-Kc + Ko)g,{-Kc + Ko))) (34) 



where Jq is decomposed according to fig. |^ for a generic SU (N) theory: 
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r{3.a) 

-'o 



Cf_ f [P + 2p' + kp) {-P + 2p' + ko){P + p + p'- kp) (P + 2p - fcp) _ 
2NJko [(f +p-A;)2-m2][(-f +p' + A;)2-m2] 

^(^-2.-) (36) 



S^) - (r^ _ 9L\ [ f iP + P + P'-ko)iP + 2p-ko) 



+ 



P + 2p' + ko){-P + p + p' + ko) 
[{-^ + p' + ky - m^] 



2(C^|-^)(^-^m) (37) 



.2 C'f . A 



= -(C'-—)- (38) 
Using the abbreviations 

f f dk^ A. 2 

L ^ J '2^ki + A^ ^^^^ 
- 1 

Cf = (40) 

2N ^ ' 

we have written the result for Pauh Villars regularization to make the cancellation of the 

liner divergent parts obvious. 

For the double Coulomb exchange graph from fig. ^ we obtain for the time component 

integral 



.(2.6) _ ^2 



[{P'o + -P^ + ^m^EkEp,] [{po + kof - P2 + Am^EkEp 
ko Hi + ko)' - El][{^ + koY - El] 



-CU^-tm) (41) 
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Collecting everything from above we have 

/o = it^ + /#"^ + /#^^ + I^'^ = -CUm, (42) 
which leads with eq. (^),(^) and immediately to the result 
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ma 



^Mabelianbox Xg^3(^/_j_l^ ^^'^^ 

The net result for the spectrum of two scalars bound by an abelian gauge field is equal 
to that of ref [Q. However, we showed which graphs contribute in a pure BS approach 
which can be used as a basis for any higher order calculation. 

We have also checked the derivative OKq/OPq contributing to hi and the X-graphs of 
fig. (§.g) with transverse gauge particles for possible contributions. Our estimates only 
yield contributions to higher order. Due to mass and wave function renormalization we 
can further assume that the graphs of fig. (^e,f) give only contributions to 0{q^ In a) as 
in the fermionic case [|12|. Possible large contributions of lighter particles to the vacuum 
polarization as depicted in (§.c) can be treated as in the fermionic case |]13| . 



In supersymmetric theories a |$|^ term is part of the lagrangian. Clearly it can be put 
in by hand into the Langragian of an ordinary qantum field theory. The contribution from 

an interaction term of the form — A/2($^T'^$)($^T"$) is easily calculated: 

■? 

AMx = -CfX^^S^o (44) 

and gives a contribution of the same form as the Darwin term (usually interpreted as a 
zitterbewegung contribution) which is suppressed by two orders in a in the scalar theory. 
There may exist a small chance that this term may be helpful for the determination of 
the supersymmety parameters of the theory contained in A. 

For an ordinary quantum field theory without a direct interaction on the tree level 
it was shown first by Rohrlich [|14| that a counter term of this form is needed for the 
scattering of two scalars (e.g the graphs of fig. § and the first of fig. Q.b with Photons in 
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Feynman gauge). It is interesting to note that in Coulomb gauge the divergencies for the 
Coulomb photons cancel and the only divergent graph is the one of fig ^.c with transverse 
photons. 

The spectrum calculated so far is common for the abelian and the nonabelian theory. 
Collecting all pieces a we have 

AM = AMl, + — 7 + ^(1 —) 771 — TT (45) 

where AMp^i originates in the contribution of j light fermions to the vacuum polarization 
and can be found in |1T3| . 

It has been pointed out first in |[T5| that in the case of a nonabelian gauge field further 



corrections may arise due to the gluon splitting vertices. The 0{a^) corrections from fig. 
(|.a, b) as well as the 0{a^) corrections from the corresponding two loop graphs are 
obviously the same as in the fermionic case. The vertex correction shown below in fig. 
(|.d) has been calculated in |jl5| for the fermionic case. 



Here we will give a calculation of the same contribution for scalar constituents. After 
performing the color trace the perturbation kernel for the second graph in fig. (^.d) reads 

rr o . 4 f (P,+p,+p'^-ko){-Po+Po+p',) ( (^)(#) \ . . 

n4.d,2 = / . — p — — — — -[pq) + — — • (4b j 

Performing the fco integration and using the scaling 

Po ^ 2m + 0(^2) 

k ^ ak 



to extract the leading contribution in a we find that 

f^m pc 



H.., = (47) 
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Adding the similar contribution from the first graph in fig. (^.d) gives 

H,,, = JJL^. (48) 

This result differs by a factor Am? from the fermionic result which is compensated by a 
corresponding difference in the wave functions to give eventually precisely the same result 
as in the fermionic case 

In view of the fact that the result depends only on the angular momentum and not on 
the spin this seems reasonable. However, we have seen in the case of the Darwin term this 
kind of reasoning sometimes fails. 

Proceeding to the graph of fig. (§.h) we observe that in contrast to the fermionic case 
the zero component integration develops Coulomb divergencies like the abelian contribu- 
tions. Since this integrations are a little bit cumbersome in dimensional regularization we 
scetch the calculation in the appendix. It turns out finally that the box graph contribution 
in fig. (|^.h) gives the same result as in the fermionc case, which was calculated recently 

((^(4./.))) = -^vr(12-7r2)(^) (50) 

The box graph with two Coulomb lines crossed, vanishes due to the color trace. Another 
box graph with the Coulomb vertices on one scalar line replaced by a seagull vertex can 
be shown to contribute to 0{a^). However, they are in principle needed to cancel the 
Coulomb singularities. 

Thus the difference in the spectrum of the scalar bound state to O(a^) compared to 
the fermionic case is entirely due to the graphs also present in the abelian theory discussed 
above. 
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4 Bound state corrections to the decay width 

Assuming that the scalar particle under consideration decays into two other particles, the 
decay width is the imaginary part of the self energy function S at the mass shell. Focusing 
on the stop quark a possible scenario could he in ^ b + Xi [0 • We shall be interested 
in terms of the order O(a^r) where F is the tree level decay width. The first part of the 
perturbation kernel due to the exact inverse propagator — — S(p^) for the bound 
state corrections to the decay width reads 

Hi = iD-^ - iDQ^ = 

14 



^ -2t{27i)^S\pi-p2)J:'im^){pl-m^){pl-m^) (51) 
In the derivation of eq. ( |^ we expanded the self energy function around the mass shell 
S(p2) = E(m2) + S'(m2)(p2 _ m^) + 0{{p^ - mY) (52) 

and we assumed that the decay width used in the zero order equation (e.g. in Dq) is given 

by 

r=-I,n5(!^. (53) 

m 



As has been first shown in |T^, the gauge dependence contained in the off shell con- 
tribution S' is cancelled by parts of the vertex correction depicted in fig. ^.f. It gives rise 
to a perturbation kernel 

H2 = Ao^(-Po +P0+ p'o), (54) 

with Aq representing the vertex correction. The color trace is already included in a. As 
in the fermionic case it is possible to derive a Ward identity which guaranties the 
cancellation of the gauge dependent terms (the 's are the SU {N) generators) 



d_ 

' dp^ 

ImAo(p = (m, 0), g = 0) = -2mImE'(m^). (55) 



A;:(p,g = 0) = -2(7T>^— E(p2 



But the detailed calculation shows differences to the fermionic case: The sum of the 
contributions from Hi and H2 vanishes to the desired order with the help of the zero 
order equation 

Im((i7i + if2)) ~0. (56) 

On the other hand, we observed above that the wave functions for decaying fermions 
and scalars were very different. While it was possible to obtain the same wave functions 
for decaying fermions and for stable ones, in the bosonic case we used wave functions 



explicitely containing the decay width (cf. sect. |2.2| ). Thus we have to reexamine the 



relativistic corrections to the energy levels. Among the contributions considered in the 
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last section only the relativistic Coulomb correction fig. (|T].a) can produce corrections to 

It is easy to see that the only difference comes from the fact that the perturbation has 
to be taken at the position of the pole ([T9|) which leads to the replacement 

al^al + z-. (57) 
m 

in eq. (^). We thus get a relativistic correction to the decay width of the bound state 

Ar,. = -^ (58) 

This has to be added to the O(a^r) term of (^9|) to yield the final result for the boundstate 
correction to the decay width: 

This result generalizes the result of 0, |1^ to the bosonic case. We can thus say that 
the effect of the bound state corrections to the decay width can be interpreted entirely as 
a time dilatation effect as was first conjectured for the fermionic theory [l^ . 



5 Conclusion 

We have presented a consistent formalism for the calculation of bound state properties 
for scalar particles interaction with an abelian or nonabelian spin one vector field. This 
is done by deriving a solvable relativistic zero order equation similar to that of Barbieri 
and Remiddi both for stabel and unstable scalars. Based on this equation a systematic 
perturbation theoy can be built which allows especially the calculation of the position of 
the bound state poles to higher orders. 

Using this approach the bound state spectrum was calculated to 0{a'^). We found 
that we had to take into account the abelian box graphs to this order. This is not the 
case in the fermionic theory. All the relativistic Coulomb corrections only reproduce the 
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p ^ term from the expansion of \/vti^ + indicating that a fully relativistic formulation is 
not really economic for the lowest orders in perturbation theory. However, the advantage 
of the presented formalism is that it is straigtformward applicable to any higher order 
calculation. We calculated also the nonabelian contributions to O(a^). Furthermore our 
approach makes possible the calculation of the bound state corrections to the decay width 
of weakly decaying scalar particles. We show that - as in the fermionic case - the inlcusion 
of a finite, constant decay width in the zero order equation simplifies the problem of the 
bound state correction to the decay width in a profound way. It is now possible to clearly 
isolate the underlying cancellation mechanism which automatically gives a gauge inde- 
pendent result which can be interpreted as time dilatation alone. We can thus gereralize 
the theorem on the bound state corrections for the decay width to the scalar case: The 
leading bound state corrections for weakly bound systems of unstable scalars (with decays 
like tfl ^ + Xi) s-re always of the form (|58D. 

It would be very interesting to observe a particle where the above mentioned predic- 
tions could be tested. Today it seems that the stop-antistop system could be a candidate. 
It will be heavy enough to allow a perturbative treatment even for the nonabelian case. 
Whether the decay width will be small enough to allow a detailed study of the spectrum 
remains open to speculation at present. But even for a quite large decay width the scalar- 
scalar potential will provide the basis for interesting threshold calculations for this case 
PPI , analogous to the ones for the top-antitop system |2T| . 
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A Zero component integrations for the nonabelian 
box graph 

The diagram ^.h leads to the energy component integrals 

dko f dto (Po + 2po - to)iPo + 2po -Qo- io)(-^o + 2po - Qo - ^o)(-^o + 2po - ko) 



/n : 



27r J 2vr [(Po/2 +po - to)' - [(-Po/2 + Po - - [(^o - koY - {t - k^] 

(60) 



This integral is power counting logarithmic divergent, but it turns out that the first 
integration is finite which leads to a linear divergent second integration. After scaling eq. 
(|60|) reduces to 

dkp r dto (2m - to)(2m + kp) 

2n J 2n (t, _ ^e) (ko + it) [(to - ^o)^ - {t- k)^ + it] 
To make this integral accessible for the methods of dimensional regularisation we use the 
following trick 

^^ = lim/^ . (62) 

2-K to ± it M-*o J 27T to - jJ,' + it 

Performing first the to integration we have 



-^0 — -^0,1 + Io,2 

_ . dko (2m + ko) {ko + 
- kl-^' + it 

dko 



/0,2 - / —.to 



2mzr(2-f) ^^{xko-ii)[x{l-x)]^-' ^r(l - f ) 



hn = n ^ / dx ^ fr—\t — k\ 

(47r)f Jo r H + ^ + ^l^-f (An- ' ' 

I ' X l—x ' 

The limit — > has to be performed very carefully to obtain 

2m^ m 

-'0,1 



|t-A;|2 2|t-fc| 

m , , 

/o,2 = -^-^ 63 
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Thus we have to the desired accuracy 

It should be noted that dimensional regularization does not show up linear divergencies. 
Instead the use of the regularization ( pQ]) leads to a visible liner divergent term (of higher 
order in a) iKj{^\t — k\) which has to be cancelled by similar contributions from graphs 
where the two Coulomb gloun vertices on one or both scalar lines are double Coulomb 
vertices. 
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